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ABSTRACT 

We present a novel chemical database for gas-phase astrochemistry. Named the Kinetic Database 
for Astrochemistry (KIDA), this database consists of gas-phase reactions with rate coefficients and un- 
certainties that will be vetted to the greatest extent possible. Submissions of measured and calculated 
rate coefficients are welcome, and will be studied by experts before inclusion into the database. Besides 
providing kinetic information for the interstellar medium, KIDA is planned to contain such data for plan- 
etary atmospheres and for circumstellar envelopes. Each year, a subset of the reactions in the database 
(kida.uva) will be provided as a network for the simulation of the chemistry of dense interstellar clouds 
with temperatures between 10 K and 300 K. We also provide a code, named Nahoon, to study the time- 
dependent gas-phase chemistry of 0D and ID interstellar sources. 



Subject headings: Physical data and processes: astrochemistry — Astronomical databases: miscellaneous 
abundances, molecules 
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1. Introduction 

The interstellar medium (ISM) is composed of gas 
and dust under a wide range of physical conditions. 
The diffuse interstellar medium is characterized by 
very low gas densities (< 10 3 cm -3 ), temperatures be- 
tween 30 and 100 K, high UV photon fluxes and chem- 
ical species mainly in atomic form (except for molec- 
ular hydrogen). Due to gravitational interactions and 
turbulence, this material slowly forms denser (~ 10 4 
cm -3 ) and colder (10 K) regions through which ion- 
izing cosmic rays can penetrate but UV photons can- 
not. Over a period of about 10 5 -10 6 yr, these atoms 
and molecular hydrogen start a series of chemical reac- 
tions that synthesize molecules in the gas and on dust 
particles. Some of these molecules are large species, 
mainly in the form of unsaturated (H-poor) carbon- 
chains. Due to perturbations, these cold and dense 
regions, called dark molecular clouds or cold cores, 
collapse to begin the different stages of star and planet 
formation; these stages comprise a prestellar core, a 
protostar, a star with a protoplanetary disk, and then 
possibly a planetary system. As the physical evolution 
occurs, a chemical evolution also occurs in which the 
dominant larger molecules are no longer unsaturated 
carbon chains but more saturated organic and prebi- 
otic species quite common on the earth. At the end of 
their lives, stars re-inject their material into the ISM, 
now enriched in heavy elements. This set of events is 
the life cycle of stars and interstellar matter. 

The study of the composition and evolution of inter- 
stellar material is part of what we call astrochemistry. 
During the next ten years, many fascinating discover- 
ies are expected in this field of research thanks in part 
to a number of new observational instruments includ- 
ing the Herschel Space Observatory (HSO), the Ata- 
cama Large Millimeter/Submillimeter Array (ALMA), 
as well as the expanded Very Large Array (eVLA), and 
the Stratospheric Observatory for Infrared Astronomy 
(SOFIA). HSO and SOFIA are currently observing 
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molecular spectra in the far-infrared, a frequency range 
not accessible from the ground, whereas the ALMA in- 
terferometer, now only in its early science stage, will 
give access to very high spectral resolution and sen- 
sitivity, with broadband spectral coverage as well as 
spatial resolution never achieved before. Many new 
molecules have been and will be detected, increasing 
our knowledge of the complexity of the ISM. More- 
over, astrochemistry is important for the understand- 
ing of Sun-like star and planetary formation. Indeed, 
protostellar collapse is governed by quantities such as 
the ionization fraction (directly related to the chemical 
composition), which couples the gas to the magnetic 
field. Finally, the evolution of the chemical compo- 
sition of the gas, dust, and ice in the protoplanetary 
disk is a key process for the understanding of plane- 
tary formation, the fate of complex organics acquired 
in earlier stages, as well as the emergence of life. 

In order to understand the formation and destruc- 
tion of molecules observed in assorted sources in 
the ISM, astrochemists model the chemical evolu- 
tion through reaction networks. In these models, the 
chemistry is described by processes of various types: 
ion-neutral and neutral-neutral reactions, dissociation 
and ionization by UV photons and cosmic ray par- 
ticles, dissociative recombination, etc. The models 
compute the abundances as a function of time for up 
to four hundred gas phase species using more than 
four thousand reactions, starting from an initial com- 
position. Extensions that include isotopologues, an- 
ions, and large molecules contain many more reactions 
and molecules. For dark clouds, the initial composi- 
tion can consist of abundances measured in more dif- 
fuse regions, mainly atoms and molecular hydrogen. 
Two main chemical databases of gas-phase chemistry 
currently exist in astrochemistry for the interstellar 
medium. The better known is the udfa database (pre- 
viously called the UMIST database), created by Tom 
Millar (Queen's University of Belfast, Northern Ire- 
land), first released i n 1991 dMillar et al.lll99l[ 1 1997 



LeTeuffet al] l2000t IWoodall et all l2007h . Updated 



versions of the database have been released in 1995, 
1999, and 2006 (http://www.udfa.net/). The second 
one is the OSU database, initially developed by Prasad 
and Huntress (1980), updated by Chun Leung and Eric 
Herbst, and then revised over the last twenty years 
by Herbst and collaborators. Several versions dating 
back to 2003 are available at http://www.physics.ohio- 
state.edu/~eric/research.html. Although the two databases 
are similar, some significant differences exist because 
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of occasionally different choices of which reactions to 
include; udfa includes more studied reactions whereas 
OSU relies more heavily on unstudied radical-neutral 
reactions. In addition, the two networks differ partially 
because of 

1 . Different approaches to estimate the rate coeffi- 
cients of reactions not studied in the laboratory, 

2. Different choices of experimental rate coeffi- 
cients from uncritical compilations (the NIST 
database for instance^, and 

3. Different approximations regarding the temper- 
ature dependence of ion-polar neutral reactions. 

Until its most recent version, which is useable 
up to 800 K, the OSU database had been more rel- 
evant for cold sources (~ 10 K) whereas udfa includes 
many reactions with energy barriers for hotter condi- 
tions (> 100 K). Other accessible databases exist for 
particular cases such a s Srat efl for sulfur chemistry 
dWakelam etail |2004 1201 ll). a datab ase for shock 
models Jpineau des Forets et al.l Il993[) and on e for 
photon-dominated region^ ( Le Petit et al. 2006). The 
existence of these diverse databases, which contain 
many unstudied reactions, and often lack references, 
is a barrier to the improvement of chemical models. 
We have therefore initiated a project of constructing a 
unique and critical database for gas-phase astrochem- 
istry named the Kinetic Database for Astrochemistrj0 
(KIDA). This database will include as many of the 
important gas-phase reactions for modeling various 
astronomical environments and planetary atmospheres 
as possible. Indeed, reaction networks relevant for 
planetary atmospheres require many of the same re- 
actions as those in interstellar modeling, but over a 
different range of temperatures. 

2. General presentation of the database 

KIDA is designed in order to allow the efficient 
managing of data, user-friendly visualization of data 
and flexible downloading of lists of reactions. In ad- 
dition to being able to consult the database, users can 
submit data to KIDA via an online interface. 



2.1. Data on chemical species 

Users can find two different types of data in KIDA. 
Some data are related to chemical species and others 
are related to chemical reactions. In the first type, 
one can find InChI codes, CAS registry numbers, elec- 
tronic states, polarizabilities, dipole moments, and en- 
thalpies of formation. InChI codes^l (for International 
Chemical Identifier) are character chains defined by 
IUPAC to identify any chemical species. The InChI 
code of OCS for instance is InChI=lS/COS/c2-l-3. 
These codes are needed for clear identification, espe- 
cially in the context of interoperability of databases. 
Free software is available to draw the structure of the 
molecule from the InChI code. CAS registry numbers 
are numeric identifiers from the American Chemical 
Society. These numbers exist, however, only for chem- 
ical species used in industry. Most of these data in 
KIDA have been collected from other databases such 
as the NIST Chemistr y WebboolU CCCBDB0, and 
IWoon&Herbs1d2009b . JPEG images of the structures 
of the species, drawn with the JMOL free software!!, 
are also displayed. 

2.2. Rate coefficients 

For chemical reactions, rate coefficients are stored 
in the form of three parameters (A, B and C), which 
can be used to compute the rate coefficients using dif- 
ferent formulae depending on the considered process 
and the range of temperatures desired. Rate coeffi- 
cients define the intrinsic velocity of chemical reac- 
tions; the units in which they are expressed depends 
on the order of the reaction. The overall rate of a reac- 
tion is given by the product of the rate coefficient and 
the densities of the reactants. There are nine classes 
of reactions in KIDA, which are listed and described 
in Table Q] The first three classes can be thought 
of as unimolecular processes induced by cosmic rays 
or photons; here the rate coefficient has units of s _1 , 
and reaction rates are computed by multiplying the 
rate coefficient by the density (in cm 4 ) of the reac- 
tant. The next five classes involve bimolecular colli- 
sions; in these classes rate coefficients have units of 
cm 3 s _1 , and the rate of reaction is computed by multi- 
plying the rate coefficient by the densities (in cm -3 ) of 



'http://kinetics.nist.gov/kinetics/index.jsp 

2 http://kida.obs.u-bordeaux 1 .fr/uploads/models/Srates_ update.txt 
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the two reactants. The last type is composed of ter- 
molecular (three-body) association reactions, which 
are processes where two reactants associate to form 
a product that is stabilized by collision with a third 
species, which can be a reactant or an inert atom or 
molecule. In the low-density limit, the rate of reac- 
tion is given by a rate coefficient in units of cm 6 s 
multiplied by the product of the densities of the three 
species involved in the reaction. A more general rate 
law involves a more complex dependence on density, 
which eventual ly saturates into a bimolecular rate law 
dAnicichlll993h . Three-body association reactions are 
effective in "very dense" sources by interstellar if not 
terrestrial standards (> 10 11 cm 4 ) and are generally 
considered in the atmospheres of planets. This type of 
reaction will be included later into KIDA. 

The rate coefficients for bimolecular processes in- 
volving two heavy species or one heavy species and an 
electron are each given a certain range of temperature 
in which the parameters are applicable. By default, 
the rate coefficients are valid between 10 and 300 K. 
For updated rate coefficients, the temperature range is 
given by the results of experiments or the estimates of 
experts. Extrapolations outside this range of tempera- 
ture are not recommended unless the rate coefficient is 
predicted to be totally independent of temperature, as 
occurs for the Langevin model of exothermic ion-non- 
polar-neutral reactions. Temperature dependent rate 
coefficients, at a specific temperature T, can often be 
fit to the Arrhenius-Kooij formula: 



k K ooij(T) = A(77300) B exp(-C/r) cm 3 s" 



(1) 



where A, B, and C are parameters. Several sets of rate 
coefficients for the same reaction can exist for different 
ranges of temperature or even the same range. If this 
happens for identified key reactions, KIDA experts can 
give recommendations on which values to use (see sec- 
tiongj). 

Rate coefficients for unmeasured reactions be- 
tween ions and neutral species with a dipole moment 
are computed using the Su-Chesnavich captur e ap- 



proach as discussed i n IWoon & Herbst (2009) and 



Wak elam et al.1 (1201 Obi) . Here the rate coefficient ku 
is expressed in terms of the temperature-independent 
Langevin rate coefficient k^ using two formulas, one 
for lower and one for higher temperatures. The tem- 
perature ranges are divided by that temperature for 
which a parameter x, defined by the relation, 



A'D 



where (in cgs-esu units) yUo is the dipole moment and a 
the dipole polarizability, is equal to two. The formulae 
for ko/kh at the lower and higher temperature ranges 
are respectively 



k D /k L = 0.4767* + 0.6200 



(3) 



and 



k D /k L = (x + 0.5090)710.526 + 0.9754. (4) 



At the lowest temperatures in the region of validity for 
equation (3) (T > 10 K), ku goes as T~ 1 ^ 2 , whereas at 
the highest temperatures, ku approaches the Langevin 
value. The rate coefficients are valid in the so-called 
semiclassical region, before they level off with de - 
creasing temperature near 10 K or below (lTroelll996h . 

Rate coefficients for photo-processes by UV pho- 
tons, kphoto, are normally computed as a function of 
the visual extinction A,,: 



kphoto = A exp(-CA v ) s 



(5) 



The rate is obtained by multiplying kphoto by the 
density of the reactant. The parameter A for kphoto 
represents the unattenuated photodissociation rate 
in a reference radiation field whereas exp(-CA v ,) 
takes into account the continuum attenuation from 
the dust. Since dust particles tend to attenuate light 
at shorter wavelengths more strongly, values of C 
are larger for species that have thresholds for pho- 
todissociation/photo ionization farther into the UV 
(van Dishoeck 1988). Photoinduced rates depend on 
the external UV radiation field assumed. For example, 
the so-called average interstellar radiation field differs 
from the radiation field around a specific type of star. 
It is possible in KIDA to specify which field will be 
used. The data included in KIDA have been com puted 
for the interstellar radiation field (Drainel ll978 ). Data 
computed with other filed will be included in the fu- 
ture. 

Cosmic rays, chiefly protons and alpha particles 
with energies above 1 MeV, interact with gas-phase 
atoms and molecules in two ways. First, these high- 
energy particles can directly ionize a species. This di- 
rect process is dominant for H2, H, O, N, He and CO. 
Electrons produced in the direct process can cause sec- 
ondary ionization before being thermalized. The sec- 
ondary process is ta ken into account in most treatments 
by a f actor of 5/3 ( Dalgarno et alJfl999T: Rimmer et al 



(lakT) 1 ' 2 ' 



(2) 



2012) while the effect of ionization by alpha particles 
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is taken account of by a separate factor of 1.8. The to- 
tal ionization rate coefficient including these two fac- 
tors is computed by 



k C R,i = Ai( 



(6) 



where ( (s ) incorporates the flux of cosmic rays 
and secon dary electrons and normally refers to atomic 
hydrogen ( Cravens & Dalgarnolll978 ) . In the udfa, 
OSU, and KIDA networks, however, f refers to 
molecular hydrogen (£ff 2 ), and is approximately twice 
the value for atomic hydrogen. The rate coefficients 
for different species can be obtained by including a 
factor A,. 

In addition to secondary ionization, the electrons 
produced in direct cosmic ray ionization can electroni- 
cally excite both molecular and atomic hydrogen. The 
radiative relaxation of H2 (and H) produces UV pho- 
tons that ionize and dissociate molecules. This indirect 
process is known as the Prasad-Tarafdar mechanism 
dPrasad & Tarafdarl 1983b . ' t represents a source of U V 
photons inside dense clouds. The rate coefficients of 
dissociation or ionization that proceed via excitation 
of molecular hydrogen are given by 



A, 



n(H 2 ) 



(l-<o) (n(H) + 2n(H 2 )) 



(H2 



(7) 



with the total ionization rate for H2, to the albedo 
of the grains in the far ultraviolet, n{H) the density 
of atomic hydrogen and n(H 2 ) the density o f molec- 



ular hydrogen (see also Gredel et al. 1989b . In the 
OSU network (and subsequently in KIDA), it is as- 
sumed that to is 0.5 and that the abundance of atomic 
hydrogen is negligible compared to molecular hydro- 
gen. The rate then simplifies to kcupj = A^h 2 ■ This of 
course is only valid in dense sources. Unlike the case 
for direct cosmic ray ionization, the factor A,- can be 
quite large, often in the range 10 3 - 10 4 . The case 
of CO m ust be treated indiv idually because of self- 
shielding dGredel et al.lll987b . 

Any data provided in KIDA is associated with a 
reference. This reference can be bibliographic, an- 
other database, or a datasheet, as discussed in Sec- 
tion 12.31 All the reactions from the low-temperature 
OSU database osu-01-2009 have be en included, even 
reacti ons involving anionic species dHarada & H erbst 
2008). Many updates concerning photoinduced rates 
are currently being under take n based on results from 



van Dishoeck et al.l (12006b and lvan Hemert & van Dishoeck 



mic ray-induced ph otons have been added based on 
Chabot et al.l Q2010). Reactions necessary for primor- 
dial chemistry, collected by D. Galli, will be included 
before the end of 2012. Finally, many other additions 
have been made, in particular for neutral-neutral re- 
actions, based on new experimental and/or theoretical 
results. All the original bibliographic references listed 
in KIDA are given in an online appendix to this paper. 

2.3. Recommendations concerning rate coeffi- 
cients 

For all the rate coefficients in KIDA, a status is 
indicated. There are four possible recommendations, 
which are listed below: 

0. Not recommended: This status suggests strongly 
that the reaction not be included. 

1. Unknown: This status is mainly for reactions 
from other databases that have not been looked 
at by KIDA experts. 

2. Valid: This is a status meant for reactions that 
are associated with a bibliographic reference. 

3. Recommended: Recommended rate coefficients 
are mainly deduced from a variety of sources 
(experimental or theoretical) but can also just 
be educated guesses. The recommended sta- 
tus must be justified by a datasheet and can 
evolve with time when new data are available. 
Datasheets are pdf files stored in KIDA0 sum- 
marizing the information available on a reac- 
tion. These datasheets are inspired from the 
work done by the IUPAC subcommittee for Gas 
Kinetic Data Evaluation^ Datasheets, of which 
there are currently more than 40, are signed by 
one or more authors and are reviewed internally 
by KIDA experts. 

2.4. Uncertainties for rate coefficients 

The reliability of any reaction rate coefficient stored 
in KIDA is expressed in terms of an estimated uncer- 
tainty over the temperature range covered by the rec- 
ommendation. This uncertainty is specified by three 
fields, two of which are shown in Table 2: 



(2008). Branching ratios for the dissociation of car- 
bon clusters by both external UV radiation and cos- 



9 http://kida.obs.u-bordeaux 1 .fr/publicationsttspan_ datasheets 
http ://w w w. iupac-kinetic . ch.cam. ac . uk/ 
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1 . The type of the statistical distribution represent- 
ing the uncertainty. The preferred distribution, 
and the default, is the Lognormal (default = 
logn). This choice enforces the positivity of 
the rate coefficient , even for large uncertainty 
factors (Fn > 2) (Thompson & Steward [l99U 

m 



Hebra rd et al Other distributions are im- 



plemented to accommodate the most common 
uncertainty representations: Normal (norm), 
Uniform (unif), and Loguniform (logu) (Table 
2); 

2. The uncertainty parameter Fq. By default, it is 
the standard geometric deviation for the Log- 
normal distribution, but has no default value. 
It defines a standard "la" confidence inter- 



val around the reference value ko (San der et al 



2001 II) . with the probability of assertion (Pr) 
defined by Pr(k /F Q < k < k * F ) - 0.68. 
For other distributions than the Lognormal, Fq 
has different meanings and units, as specified in 
Table 2; 



3. The expansion parameter g (default = 0, units: 
Kelvin). It is used to parameterize a pos- 
sible temperature-dependence of the uncer- 
tainty, according to t he formula F(T) - Fq * 



exp(gll/r- 1/TqI) (Sa nder et al] 12001 ll) . In 



KID A, the reference tempe rature To is fixed at 
300 K dHebrard et al.ll2009h . 

More details on how to use, choose, and specify un- 
certainty representations in KIDA, notably for data 
providers, is given in http://kida.obs.u-bordeauxl.fr/ 
help/helpUncert.html. 

2.5. Member Area, Submitting Data to KIDA, 
and Downloading Lists of Reactions 

In order to be able to interact effectively with KIDA, 
users need to first register by providing an email ad- 
dress and a password. They will then have access 
to a member area where they can submit rate coef- 
ficients to KIDA experts and see the history of their 
downloads. To improve the updating of KIDA, we 
have designed an online interface to upload data. Any 
user registered with KIDA can submit rate coefficients 
(with a temperature range of validity, an uncertainty 
and a reference). These data will then be sent to 
the "Super-Expert" of the type of reaction, who will 
then consult his/her experts to see if the data have 



been correctly included or not. Lists of KIDA super- 
experts and experts with terms until 2013 are given in 
http://kida.obs.u-bordeauxl.fr/contact 

It is also possible to download lists of reactions 
from KIDA using an online forrrFl. In this form, the 
user can select 1) the elements and/or species he is in- 
terested in, 2) the types of reactions and 3) the range of 
temperature. This last choice can be used, for instance, 
as a filter to remove some of the high temperature re- 
actions if one is interested in modeling the cold inter- 
stellar medium. If no recommended rate coefficients 
exist in the full range of temperature asked by the user, 
several rate coefficients for the same reaction may be 
available in the downloaded file. These reactions are 
included in the readme file downloaded at the same 
time, and the user has to deal with them. If rate coeffi- 
cients are defined over complementary ranges of tem- 
perature, the user may want to keep them and take into 
account this information in the model (see section lXTl i. 
Temperature ranges are indicated in the network along 
with the status of the KIDA recommendation (see sec- 
tion |33]l. In any case, we do not recommend the ex- 
trapolation of rate coefficients outside this range. If 
this procedure must be undertaken, we recommend in- 
stead using the results at the relevant extremum of tem- 
perature or consulting an expert or super expert from 
the KIDA lists. 

For each reaction, KIDA stores up to three reac- 
tants, up to five products, the values of the kinetic pa- 
rameters A, B and C, the uncertainty parameters Fq, 
g and the type of uncertainty, the temperatures T m i n 
and T max , which define the range of temperatures over 
which the rate coefficient is defined, the name of the 
formula used to compute the rate coefficient or the 
external UV field used to compute the photoinduced 
rates, the number of the reaction in the network, an- 
other number that represents a particular choice of A, 
B and C when more than one possibility has been 
included, and a numerical recommendation given by 
KIDA experts and discussed in the next section. When 
there is no relevant temperature range, -9999 and 9999 
default values are used for the minimum and maximum 
temperatures. 



http://kida.obs. u-bordeaux 1 .fr/astrophysicist 
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3. New Model for the Cold Dense Interstellar 
Medium 

3.1. The Nahoon Public Code 

The numerical code Nahoon is publicly available on 
KIDA to compute the gas-phase chemistry for astro- 
physical objects (http://kida.obs.u-bordeauxl.fr/models/). 
The program is written in fortran 90 and uses the DL- 
SODES (double precision) solver from the ODEPACK 
package (http://www.netlib.org/odepack/opkd-sum) to 
solve the coupled stiff differential equations. The 
solver computes the chemical evolution of gas-phase 
species at a fixed temperature and density (commonly 
referred to as zero-dimension, or OD) but can be used 
in ID if a grid of temperature, density, and visual 
extinction is provided. Nahoon is formatted to read 
chemical networks downloaded from KIDA. It con- 
tains a test to check the temperature range of the va- 
lidity of the rate coefficients and avoid extrapolations 
outside this range. A test is also included to check 
for duplication of chemical reactions, defined over 
complementary ranges of temperature. In that case, 
Nahoon chooses the recommended rate coefficient. 
Some reactions, although important for the interstellar 
medium, are not included in KIDA such as H2 forma- 
tion on grains or recombination of cations with neg- 
atively charged grains. These reactions are, however, 
provided in a set of user instructions associated with 
Nahoon. This code represents a new version of the 
one t hat was put online in 2007 (Wa kelam & Herbst 
2008). Grains, both neutral and negatively charged, 
and electrons are considered as chemical species and 
their concentrations are computed at the same time as 
those of the other species. 

3.2. Gas-phase network for cold dense molecular 
clouds 

KIDA is a non-selective database. It may not be 
easy for all users to select the most appropriate list 
of reactions. For this reason, we will release every 
year a standard list of reactions for the dense in- 
terstellar medium between 10 and 300 K. This net- 
work will be available on the same page as Nahoon 
(http://kida.obs.u-bordeauxl.fr/models). To treat the 



formation of molecular hydrogen on granular sur- 
faces, the OSU database uses the approximation from 
Hollenbach & Salpeter Jl97lh . which neglects any de- 
tails of the process, and contains the assumption that 
the formation rate of H2 is equal to 1 /2 the colli - 
sion rate of H atoms with grains (see also iLeung et"al 



1984). This approximation can induce some numerical 
problems during the computation because the elemen- 
tal abundance of hydrogen is not rigorously conserved 
by the rather complex procedure used by the authors 
of the code. Rather than improving the procedure, 
we use a simplified version o f the more detailed ap - 
proach in gas-grain models dHasegawa et alJ [l992). 
In particular, we have divided the H2 formation into 
two reactions and included another species JH, which 
refers to hydrogen atoms on grains. The first reaction 
is the sticking of atomic hydrogen onto grains (neutral 
or negatively charged): H — > JH. 

The rate of this reaction (in cgs units of cirT 3 s ) 

is 



R 



H,acc = 



8/tr 



-n d n(H) 



(8) 



with a the radius of grains (in cm), k the Boltzmann 
constant, m# the atomic mass of hydrogen, n d the total 
density of grains, T(K) the temperature, and n(H) the 
density of atomic hydrogen. For 0. 1 /mi radius grains, 
the rate is 7.92 x 10~ 5 y/T /3Q0n d n{H). 

The second reaction is the formation of H2 on 
grains: JH + JH — > H2. We assume here that H2 formed 
on grains is directly released in the gas phase. Since 
desorption of JH is not included, the process has an ef- 
ficiency of 100%. The rate of this reaction is given by 
the expression 



R 



ypexp (-E b /kT) 2 

H 2 ,form - — n(JH) 

N s n d 



(9) 



with N s the number of binding sites on a grain, vq the 
characteristic or trial vibrational frequency, and Et the 
potential energy barrier between two adjacent surface 
potential wells. 

A useful approximate formu la is vq = ^2n s EukBl(nmH) 
( Tiel ens & Alla mandola 1 19871) . where n s is the sur- 
face density of sites on a grain in cm 2 , and Ed 
the adsorption energy of H in er g. The read er 



can fi nd a detailed explanation in lHasegawa et al. 



(1992) and ISemenov et al. I d2010h . Assuming grains 
of 0.1 //m radius, the num ber of sites per grain 
of 10 6 dHasegawaetailll 992b . E D /k = 450 K, and 
E b lk = 225 K, n s is 7.96 X 10 14 cm" 2 , v is 4.32 x 
10 12 s _1 and the rate of H2 formation is 8.64 X 



10 b ex p (-225 /T)n( J H) z /n d cm 3 s _1 . See lGarrod & Herbsl 



(2006) for a discussion concerning the values of 
Ed and Et,, which depend upon the surface. The 
value of n d should obey the constraint that na 1 n d m 
1 - 2 X 10~ 21 « H , so that here n^5x 10~ 12 « H , where 
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a is the granular radius and «h is the total hydrogen 
nuclear density. 

The second class of processes involves electron at- 
tachment to grains and the recombination of cations 
with negatively charged grains. Rate coefficients for 
such processes can be obtained in the ballistic approx- 
imation from the following equation: 



knntt - na 



%kT 



(10) 



where m,- is the mass (g) of the impacting species (elec- 
trons or ions). We make the approximation that the 
mass of the impactor is negligible compared to the 
grain mass, so that the reduced mass is then approx- 
imately the mass of the lighter collider. For collisions 
between two charged species (cations and negatively 
charged grains), this rate can be multiplied by (1 + ^) 
to account for the Coulomb attraction, with e the ele- 
mental electronic charge a nd Z the number of nega - 
tive charges on the grain dWeisheit & Uphamlll978b . 
It is typically assumed, however, that, in dense clouds, 
grains can be neutral or singly negatively charged only, 
reflecting the possibility that the trajectory of an elec- 
tron neari ng a negatively-cha rged grain is repulsive. 
Following [Leung et al.l (119841). we neglect any surface 
electronic charge, based on the assumption that the 
negative charge/charges on a grain remain localized so 
the additional factor does not play an important role 
in the dynamics of the recombination with cations. If 
the negative charges on a grain were delocalized, then 
the Coulomb interaction should be included. We have 
listed in Table [3] the minimum set of reactions with 
grains that should be included in a network for dense 
clouds. The rate coefficients (cm 3 s _1 ) for these pro- 
cesses are computed from the A and B listed in the 
table using the equation k = A(T /300) B , where A is 
obtained for a grain radius of 0. 1 fim and B = 0.5. For 
grains of different radii, the ballistic rate coefficients 
simply scale as the square of the radius. The referee 
has noted correctly that reactions involving gaseous 
ions and neutral grains are not included in the KIDA 
network. They are not included because the only reac- 
tive gas-grain processes in the network involve charge 
neutralization of negatively-charged grains. 

In the interstellar medium, molecules can be dis- 
sociated by absorption of UV photons via electronic 
continua or through lines. The photodissociation of 
the dominant species H2 and CO is dominated by line 
absorption (although the mechanisms are different) so 
that H2 and CO self-shielding is an important process. 



To avoid a detailed treatment of radiative transfer , 
we adopted the approximation from lLee et al. (1996), 
which gives H2 and CO photodissociation rates as a 
function of A v , as well as the H2 and CO column den- 
sities. To include self-shielding, the user of Nahoon 
has to set the conditions at the border of the cloud, 
including the abundances of H2 and CO even when 
working in 0D. In ID, a plane parallel geometry is as- 
sumed and the UV photons enter the cloud from one 
direction. Photodissociation rates at a spatial point are 
then computed using the column densities of H2 and 
CO according to the computed abundances and phys- 
ical structure in front of this point. The visual extinc- 
tion is then computed using the H2 column density and 
the standard ratio between this column density and vi- 
sual extinction. This geometry is not correct for clouds 
where UV photons can penetrate from both sides of the 
cloud. Note t hat new sh i elding functions have been 
published bv lVisser et al.l (12009b and will be included 
in Nahoon in the future. 

In addition to these processes, which are down- 
loaded with the new Nahoon, we have downloaded 
our first annual network of gas-phase chemical reac- 
tions necessary to compute the chemistry in the inter- 
stellar medium where the temperature is between 10 
and 300 K and where surface chemistry is not critical. 
Approximations for the depth-depende nt photodisso- 
ciatio n rates are valid for A v above 5 (Ivan Dishoeckl 
1988). The full network, designated kida.uva.201 1, 



and the associated list of species are available at 
http://kida.obs.u-bordeauxl.fr/models This list con- 



tains 6090 reactions involving 474 species, including 
neutral and negatively charged grains, as well as elec- 
trons. Elements that are included are H, He, C, N, 
O, Si, S, Fe, Na, Mg, CI P and F. Anionic chem- 
istry from lHarada & Herbstl d2008l) has also been in- 
cluded. For temperatures in the range 300 K - 800 
K, modelers can use the high-temperature OSU net- 
work osu-09-20 1 0-ht, found and discussed on the URL 
http://www.physics.ohio-state.edu/~eric/research.html 

3.3. Using uncertainties for rate coefficients in 
models 

Using uncertainty information in chemical model- 
ing is vital, notably for the extreme environments tar- 
geted by KIDA, where many reaction rate coefficients 
are known with little accuracy, either estimated or ex- 
trapolated. Uncertainty management has two goals: 

1 . Uncertainty Propagation: to estimate the preci- 
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sion of the model outputs; and 

2. Sensitivity Analysis: to identify key reactions, 
i.e. those contributing notably to the uncer- 
tainty of model outputs and for which better ex- 
perimental or theoretic al estimations are needed 
(Dob rijevic et alJl2010h . 



Improvement of data in KIDA is based on an iterative 
strategy involv ing uncertainty propag ation and sensi- 
tivity analysis ( Wake lam et alj|2010t 



Monte Carlo Uncertainty Propagation (MCUP). 

The simplest way to implemen t these methods is 
throu g h Monte Carlo sampling ([Thompson & Stewart 
19911: iDobriievic et all kboa IWakelametalJ 120051 
Carrasco et al.ll2007t IWakelam et al.ll2010bl) . Random 
draws for a rate coefficient at any temperature are gen- 
erated using the formulae in Table [4] For MCUP, the 
code is run for N random draws of the M rate constants 
of the chemical scheme {k^; i = 1,M; j = l,N] (all 
parameters are perturbed simultaneously), and the N 
sets of outputs are stored for further statistical analysis 
(mean value, uncertainty factor, input/output correla- 
tion...). 

Sensitivity Analysis. A convenient way to perform 
sensitivity analysis in the MCUP setup is to calculate 
the correlation coefficients between the inputs (£,) and 
outputs of the model. Large absolute values of corre- 
lation coefficie nts reveal strong influences of inputs 
on ou tputs dWakelam et al. 2010b : Dobrrjevic et al. 
l2010h . Variations of this method include using rank 
correlation co efficients, or the l ogarithm of input s 
and/or outputs dHelton et al.l200dlSaltelli et all2 004). 



Monte Carlo simulations have been implemented in 
the Nahoon code bundle, and an IDL procedure has 
been provided to compute the standard deviation for 
the predicted abundances of species due to uncertain- 
ties in reaction rate coefficients. 

3.4. Limitations 

The interstellar chemical model defined by both 
the numerical program Nahoon and the osu.uva.2011 
network takes into account a limited number of pro- 
cesses, which makes its use valid only over a limited 
range of temperature, density, and visual extinction. 
This model is meant to compute the gas-phase chem- 
istry in regions where the visual extinction is larger 
than m 3 — 4. For smaller A v , detailed PDR codes, 



such as the Meudon code dLe Petit et al.1 120091) . are 
more appropriate. Also, this model does not take into 
account positively charged grains, which play a ma- 
jor role at low extinction. The temperature limit is 
set by the network to the range 10 to 300 K, as al- 
ready discussed. For higher temperatures, one should 
adapt an expanded gas-phase network developed by 
Harada et al. (2010), which can be utilized up to gas 



temperatures of 800 K. The accretion time scale onto 
grains of radius 0.1;U is roug hly 10 9 /n H (cm- 3 ) yr. At 
times larger than this, the gas-phase can be signifi- 
cantly depleted of heavy species if the temperature lies 
below the the evaporation temperature of easily va- 
porized species such as CO and CH4. At lower tem- 
peratures, non-thermal desorption mechanisms, such 
as photodesorption, can still maintain a residual gas 
phase. Models includin g interactions with grains and 
grain surface processes (IGarrod et al.ll2.007l) are, how- 
ever, based on many uncertain processes. In addition 
their use and the analysis of the results are more diffi- 
cult. For this reason, pure gas-phase models can still 
be useful even for dense cloud conditions but users 
have to be careful when they compare their results with 
observations. 

The computation of PL formation on grains and 
the recombination of negative grains with cations, pre- 
sented here, depend on the size of the grains. We have 
computed the rate coefficients for these processes for a 
single grain size of 0.1 /urn, which with a typical grain 
density of 2.5 gm/cm 3 , can explain observed extinc- 
tions with a gas-to-dust mass ratio of about 0.01, a rea- 
sonable value con sidering heavy elemental depletions 
in diffuse clouds (lTielensll2005l) . The effects of using 
a detailed grain size distribution on the rates of surface 
processes and on the chemic al results for a full gas - 
grain model can be found in lAcharvva et al.1 (1201 lb : 
these authors included grains large enough to have a 
constant rather than stochastic temperature and used 
a var i ety of size distribu tions dWeingartne r & Draind 



200 ll: iMathis et al.lll977l) . They found little effect on 



the overall chemical results for a cold core with the 
OSU gas-grain network. The relationship between 
granular size and rates of assorted processes is a com- 
plex one, because the rates can depend upon a num- 
ber of parameters in addition to surface area, such as 
whether the system lies in the accretion or diffusion 
limit and, for cold dense cores, the depth of the ice 
mantle. 

The temperature dependence of the rate of PI2 
formation on grains is a rather complex issue, and 
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depends strongly on how rough the dust surface is 
and the availability of chemisorption sites, both of 
which can allow the process to continue up to much 
higher temperatures than the standardly-assumed 
Langmuir-Hinshelwood mechanis m based on ph- 
ysisorption on smooth surfa ces dCuppen & Herbst 
2005; Cazaux & Tielensl 2010). The net result is that 
for systems with both physisorption and chemisorption 
sites, the efficiency of H2 formation can remain high 
at all relevant interstellar dust temperatures, although 
this result is not guaranteed. For small particles, in 
which the temperature is stochastic, the situation is 
more complex bu t has also been stu died via Monte 
Carlo techniques (Cup pen et al. 2006). For dust par- 
ticles with well-defined temperatures up to 50 K, the 
physisorption mechanism assumed here should be ad- 
equate for rough particles. 

Finally, despite their possible impact on the inter- 
stellar chemistry o f cold dense cores, as shown by 
IWakelam & Herbsj (2008), Polycyclic Aromatic Hy- 
drocarbons (PAHs) are not included in this model be- 
cause of the strong uncertainties that still remain (1) 
on the size and abundance of these species in dense 
regions and (2) on their reactions with other species. 

3.5. Comparison with other networks 

To compare the abundances computed with the 
kida.uva.201 1 network and other public networks, we 
have considered the case of dark molecular clouds us- 
ing Nahoon in 0D. The temperature is set to 10 K, 
the proton density to 2 x 10 4 cm" 3 , the total H2 
cosmic ray ionization to 1.3 x 10" 17 s" 1 and the vi- 
sual extinction to 30 . Initia l abundances are taken 
from [Wakelam et al.l (I2010al) . The results for se- 
lected molecular abundances, as computed with the 
kida.uva.201 1 list of reactions, are displayed in Fig.Q] 
with the results of the uncertainty propagation of 
rate coefficients. For comparison, abundances com- 
puted using the udfa06 network (http://www.udfa.net/) 
and osu-01-2009 network (http://www.physics.ohio- 
state.edu/~eric/research.html) networks are also in- 
cluded. Contrary to the publicly available udfa06 net- 
work, the OSU network used also contains anions. 

The particular molecular abundances plotted have 
been selected because they show some variation in 
their abundances with the database used. The OSU and 
kida.uva networks lead to similar results since KIDA 
is mainly based on the OSU compilation although a 
number of updates have been performed. Abundances 



computed with udfa06 on the contrary can be different, 
as shown especially for CS and the cations in Fig. Q] 
Model predictions done with the different networks are 
however within the 2cr error bars computed from the 
uncertainty propagation of rate coefficients, which are 
typically an order-of-magnitude in both directions. 

4. Conclusions 

We have presented in this article a new chem- 
ical database for gas-phase astrochemistry named 
KIDA (for Kinetic Database for Astrochemistry). 
KIDA can become a useful tool for physical chemists 
to circulate their relevant work for astrochemistry, 
especially their most up-to-date data. In addition 
to kinetic data for the chemistry in the interstellar 
medium, KIDA will soon contain kinetic data for 
planetary atmospheres (other than the Earth). Ki- 
netic data for carbon-rich and oxygen-rich circum- 
stellar shells are also planned for the future. State- 
to-state reactions, excited species, isotopologues and 
ortho and para forms will be considered for addition 
in KIDA. The quality of the data present in KIDA 
will be improved over the years by KIDA experts and 
will be preserved thanks to the detailed information 
stored in KIDA. The visibility of KIDA will be im- 
proved by its interoperability with the Virtual Atomic 
and Molecular DataCenter (Dubernet et al. 



2010, 



http://www.vamdc.eu/) and the Europlanet Research 
Infrastructure (http://www.europlanet-eu.org/outreach/). 
In addition to the rate coefficients, we are providing the 
community with a public fortran code, able to read net- 
works downloaded from KIDA, called Nahoon, which 
can calculate the evolution of the gas-phase composi- 
tion for 0D and ID sources in the ISM and an extensive 
list of chemical reactions for the dense ISM with tem- 
peratures between 10 and 300 K. Neither KIDA nor 
Nahoon treats the interaction of the gas-phase with 
interstellar grains, except for the formation of molec- 
ular hydrogen and recombination of the major cations 
with negatively charged grains. Omission of these 
processes is not realistic in the context of interstellar 
chemistry but the numerical and chemical complexity, 
and additional uncertainties brought by the treatment 
of surface chemistry, make pure gas-phase models still 
needed. In addition, gas-grain models also rely on the 
quality of gas-phase rate coefficients. 
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A. Appendix 



A.l. KIDA references 



Orig inal references for kin etic data in KIDA are: lAdams & Smithl d 197611 1 984h:lAmand d 199dh: Anicich & Huntress 



(1986); A nicich et al.1 dl976l): 



Antioov et al. 



d2009h: lAtkinson et all d2004: IChenel et all d2010h: iBarckholtz et al 
(|2001 | ) : |Barlow & Punnl(ll987l):lBaulch et alJd2q05h:lBecker et al.ldl989l.l200d):lBedianian et al]dl999b:lBergeat et al 



(2009): lBerkowitz et al.ldl986b:lBerteloite et al.ld2010blah:lBettens et al.ldl999h:lBettens & Herbsddl995h:lBohme et aj. 



1988); Br own et all dl983h: iBrownsword et all dl997lj 1996b; ICanosa et al.1 dl997h : ICartv et all i20(m. IChabot et all 
(120101): ICostes et all d2009l): baueev et al I d2005h : bavidson & Hansonl dl99dh: beevamulla & Husainl d200< 
|Pefrees& McLean! (U989h : IPerkatch et al.1 dl999h; borthe et al.1 dl99lh: [Puff & Sharmal dl996[) : lEhlerding et all 
d2004l): lEichelberger et ail d2007l): iFehsenfeld et all dl975l): jGamallo et all d2006l): iGeppert et all d2004bllal l2000h: 
iGrebe & Homann|jl982h;lGredel et al.ldl989h:lGu et al.ld2Q06l):ISun & Lild2004:lHarada & Herbsd (bq08l):lHerbst et al 
dl989ah: iHerbsd dl985l Il987l): iHerbst et alJ dl989bl 120001): Ijensen et all d2000h: iKaiser et al.1 d2002l): iKalhori et all 
( 2002|) ;| Lee et al.l d 1 978b : iLilenfeld & RichardsoDldl977l):|Liaue et alJd2q09l):lLoison & Bergeatld2009h ; lMarciuette et all 
( 198a) : lMartinet et al.ld2004:lMartinez et alJd2qq8[l2010|):lMaver et alJdl967h:lMcGowan et alJdl979l):|Mebel & Raise 
( 2002) : iMessing et al.l*(fl98Qh : iMillar & Herbsd dl990h: Montaigne et alJ d2005b: iNeufeld et alj d2005h : bzeki et al. ~ 
( 20o4: IParisel etalJ dl996b : IPrasad & Huntress! dl980): iHuo etalJ d201 lb: ISablier & Fuiiil d2002h: ISander et al 
32001 lb:IShin & Beauchamd dl986l): ISims et all Jl994:lsmith & Adamsl dl98lb : ISmithl dl989j) : ; ISmith et all d2001 ) 

;etalJ d2010b : jTalbi & Smithl d2009b: iThorneetalJ (Il983l \19s4i: lTunaetal.1 J 2OO8) 
van PishoecM (11988b: IVan Harrevelt et alJ d2002b: Ivan Hemert & van PishoeckJ 02008 ) 
Vazquez et al.1 (l2009h: IViggiano et al.l d2005b: IVinckieil dl979b : IVroom & de Heed dl969b : IWakelam et alJ 12009); 



Wennber g et al.l dl994 : lwoon & Herbsd d2009): ZabarnicketaD 0989). 
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Table 1 

Types of reactions defined in KIDA 



Type of reaction 


Abbreviation 


Description 


Direct cosmic ray processes 


CR 


Dissociation or ionization of species due to direct collision with cosmic ray particles. 


Photo-processes induced by cosmic rays (secondary photons) 


CRP 


Dissociation or ionization of species due to UV photons emitted following H 2 excitation. 


Photo-processes 


Phot 


Dissociation or ionization of neutral species by UV photons with a standard interstellar UV field. 


Bimolecular reactions 


Bimo 


Neutral-neutral (A + B -> C + D), ion-neutral (A + + B -> C + + D, A" + B -» C" + D) a , 
anion-cation (A + + B~ — » C + D) reactions and dissociative neutral attachment (A + B — > AB + + e~) 


Charge exchange reactions 


CE 


A + + B -» A + B + and A + + B~ -> A + B 


Radiative associations 


RA 


Association reactions between two species (neutral or ionized) stabilized by the emission 
of a photon (A + B — > AB + photon or A + + B — > AB + + photon). 


Associative detachment 


AD 


Association of a neutral species and an anion, resulting in the ejection of the extra electron (A - + B — > AB + e~). 


Electronic recombination and attachment 


ER 


Dissociative neutral recombination (AB + + e~ — > A + B) or the emission of a photon 
(AB + + e~ — > AB + photon) or the attachment of an electron (A + e~ — » A~ + photon). 


Third-body assisted association 


3-body 


Association reactions between two species, stabilized by collision with a third body (A + B + C -» AB + C). 



"Charge exchange reactions have been defined as a separate type due to the particular microphysics behind this process. 



Distribution Formula 

Lognormal k(T) = exp(rn&o(r) + lnF(r) * N(0, 1)) 

Normal k(T) = lc Q (T) + F(T) * N(0, 1) 

Uniform k(T) = k (T) + F(T) * 2 * (U(0, 1) - 0.5) 

Loguniform k(T) = exp (In k Q (T) + In F(T) * 2 * (t/(0, 1) - 0.5)) 

Table 4: Generating random samples for MCUP. 
U(0, 1) is a standard uniform random number gen- 
erator (between and 1); A^(0, 1) is a standard nor- 
mal/gaussian random number generator (centered at 
0; variance 1). In systems with a temperature gradi- 
ent, and where multiple uncertain rate coefficients are 
managed simultaneously, care has to be taken that the 
same random number is used for the whole tempera- 
ture range of a single reaction. 
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Table 2: The types of distributions and uncertainty factors used in KIDA. 



Type 


Distribution 


Fq meaning 


Fq unit 


norm 


Normal 


stdev 1 : Pi 2 (k Q - F < k < k + F ) ~ 68 % 


as rate coefficient 


unif 


Uniform 


half range: Pv 2 (k -F <k<k + F Q ) = 100 % 


as rate coefficient 


logn 


Lognormal 


geometric stdev 1 : Pi 2 (k Q /Fo < k < k Q * Fq) ~ 68 % 


no units 


logu 


Loguniform 


geometric half range: Pr 2 (fco/^o < k < ko* Fq) = 100 % 


no units 



1 stdev = standard deviation 

2 Pr = Probability of assertion 



Table 3: Reactions involving grains 



Reaction 




A (cm 3 s 


- 1 ) 


B 


C + + GRAIN" 


-> C + GRAINO 


2.3 x 10" 


-5 


0.5 


Fe + + GRAIN" 


Fe + GRAINO 


1.1 x 10" 


-5 


0.5 


H + + GRAIN" 


— > H + GRAINO 


7.9 x 10~ 


-5 


0.5 


He + + GRAIN 


^ He + GRAINO 


3.9 x 10- 


-5 


0.5 


Mg + + GRAIN 


- ^ Mg + GRAINO 


1.6 x 10" 


-5 


0.5 


N + + GRAIN" 


-> N + GRAINO 


2.1 x 10- 


-5 


0.5 


Na + + GRAIN 


- -> Na + GRAINO 


1.6 x 10- 


-5 


0.5 


+ + GRAIN" 


O + GRAINO 


2.0 x 10- 


-5 


0.5 


S + + GRAIN" 


-> S + GRAINO 


1.4 x 10" 


-5 


0.5 


Si + + GRAIN 


Si + GRAINO 


1.5 x 10" 


-5 


0.5 


H+ + GRAIN 


— > H 2 + H + GRAINO 


4.6 x 10" 


-5 


0.5 


HCO + + GRAIN" -> H + CO + GRAINO 


1.5 x 10" 


-5 


0.5 


e- + GRAINO 


-> GRAIN 


3.4 x 10" 


-3 


0.5 
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Fig. 1. — Abundances as a function of time for selected species computed for dense cloud conditions (see section [331 ) 
using three different chemical networks: udfa06 (blue dashed lines), osu-01-2009 (red dashed-dotted lines) and 
kida.uva.201 1 (black lines). Dotted curves represent 2<x error bars on the abundances computed with the kida.uva.201 1 
network using uncertainty propagation in the reaction rate coefficients. 
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